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1  Summary 

Cloud  features  have  an  important  impact  on  jglobal  Navy 
activities,  ranging  from  the  effects  of  clouds  on  surface  ship 
and  aircraft  operations  to  the  limitation  and  enhancement  of 
surveillance  activities.  We  have  been  developing  an  automated 
system  for  the  recgnition  of  operationally  important  cloud  types 
using  data  obtained  from  the  Defense  Meteorological  Satellite 
Program  (DMSP)  satellite  system.  The  satellite  is  deployed  in  a 
sun-synchronous  morning  orbit,  and  carries  a  variety  •  - 

of  sensors,  including  the  Optical  Line  Scanner (OLS) .  This  paper 
describes  ongoing  work  in  the  use  of  physical  and  textural 
measures  for  the  automatic  recognition  of  cloud  classes.  We  have 
found  the  physical  measures  to  be  relatively  computer  intensive, 
and  dependent  on  visible  and  IR  satellite  data.  The  addition  of 
textural  measures  shows  promise  for  enhancing  the  classification 
capability  by  ucir.g  only  a  aingie  channel  and  by  reducing 
computer  processing  time. 

2  DMSP  Data  analysis 

Satellite  image  data  from  the  DMSP  Optical  Lines  Scanner  are 
received  at  FNOC  within  two  hours  of  observation.  The  OLS  data 
are  measured  in  two  spectral  bands,  visible  and  thermal  infrared, 
with  constant  ground  resolution  across  the  satellite  swath.  The 
data  are  immediately  transformed  to  the  Quadrilateralized 
Spherical  Cube  coordinate  system  (Chan  and  O'Neill,  1978) .  This 
coordinate  system  is  particularly  suited  to  objective  image 
analysis,  since  the  coordinate  system  is  area  preserving, 
providing  a  constant  ground  resolution  anywhere  on  the  globe. 

The  satellite  data  are  stored  in  64  by  64  pixel  interpolation 
blocks  (IB).  overlaps  front  different  orbits  are  permitted  amena 
IB's,  but  within  each  IB  data  are  guaranteed  to  be  obtained  from 
the  same  orbit. 

The  coordinate  system  is  ideal  for  multispectral  analysis  since 
the  data  are  located  on  a  equal  area  projection,  with  the  same 
resolution  for  visible  and  IR  channels 


3  Image  feature  derivation 


The  digital  imagery  are  analyzed  within  each  interpolation  block. 
The  earlier  physical  image  analysis  (  Garand,  1988,  Goroch,  1988) 
have  been  found  to  be  computer  intensive,  and  inefficient  on  the 
Cyber  175  computer  system  in  use  at  FNOC.  The  analysis  was 
expanded  to  include  textural  measures  (  Welch,  1986)  in  an  effort 
to  provide  a  graded  system  of  cloud  classification.  The  current 
objective  is  to  analyze  the  image  with  progressively  more  complex 
classification  measures,  until  a  required  confidence  level  is 
attained . 

4  Physical  features 

The  physical  features  of  the  visible  and  infrared  image  have  been 
described  by  Garandand  Weinman (  1986)and  Ebert  (1988)  and  are 
summarized  in  table  1.  The  features  are  combined  into  an  image 
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vector  which  is  normalized  and  used  in  later  feature  discrimina¬ 
tion.  The  physical  meaning  of  the  feature  vector  is  determined 
by  its  use  in  the  function  discriminating  among  cloud  types.  In 
the  current  use  of  a  Gaussian  multivariate  type  of  discriminant 
function,  the  elements  of  the  feature  vector  of  a  category  pro¬ 
vide  the  following  information 

Element  negative  -  feature  reduces  class  membership 

Element  zero  -  class  membership  distributed  around  mean 

Element  positive  -  feature  increases  class  membership 

The  Garand  cloud  classes  were  obtained  from  a  training  set  of 
GOES  data  during  the  winter  of  1984  in  the  northwest  Atlantic. 

The  classifications  are  shown  in  Table  2,  together  with  objective 
characteristics  of  each  cloud  type.  Figure  1  shows  the  charac¬ 
teristic  feature  vectors  associated  with  low  layered  cloud  types. 
Note,  "for  example,  the  improtance  of  the  MC,  multilayer,  element 
in  differentiating  the  altocumulus  cloud  type  from  the  two  multi¬ 
layered  cloud  types,  Cumulus  with  Altocumulus  and  Stratocumulus 
with  Altocumulus. 


Table  1 

Physical  characteristic  vectors 


Crosscorrelation  ( ST  .  ) 

(  0, 

l)  | 

Spectral  Intensity  between  20-40  km 

!  0 , 

i)  1 

!  SE  ) 

ll 

ne;-'.  cloud  albedo  (AL) 

S  0, 

l) 

Cloud  top  height(HT) 

0  - 

14 

km 

Background  connectivity  ( BC ) 

<  0, 

1) 

IR  low  cloud  fraction  (LO) 

<  0, 

1) 

IR  middle  cloud  fraction  (MI) 

(  0, 

1) 

Multilayer  index  (ML) 

(  0, 

1) 

Total  cloud  fraction  (CF) 

(  0, 

1) 

Number  of  clouds  ( NC ) 

Fraction  of  cirrus  clouds  ( LR ) 

(  0, 

1) 

Cloud  Connectivity  (CC) 

(  0, 

1) 

High  Cloud  fraction  (HI) 

■  'a 

5  Textural  features 

To  supplement  the  physical  features,  we  have  added  the  textural 
features  used  by  Welch!  1986)  and  others  for  class  discrimina¬ 
tion.  The  textural  features  are  shown  in  Table  3.  The  textural 
features  provide  a  rapid  textural  analysis  which  it  was  consid¬ 
ered  would  add  to  useful  measures,  which  would  not  use  every  much 
computer  time. 
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Table  2  -  Physical  class  characteristics 
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Cloud  Feature 

Figure  1,  Low  layer  cloud  feature  vectors.  Altocumulus ( square )  , 
Cumulus  with  Altocumulus  (plus),  Stratoeumulus  with  Altocumulus 
( d iamond ) 


Table  3  Textural  Feature  Analysis 
Grey  Shade  Co-Occurrence  Calculation  Functions 
nv{i.p).v(i.P-d)) 

/(V(L.p).V(l-d.p)) 


where  f  =  sum,  difference,  Absolute  value  of  difference 
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Grey  Shade  Co_occurence  Texture  Measures 
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The  textural  features  were  calculated  by  determining  the  frequen¬ 
cy  of  occurrences  of  differences,  sums  and  absolute  differences 
in  pixel  values  between  pixels  separated 'by  d  pixels.  The  current 
analysis  was  conducted  for  separations  of  d-0,  adjacent  cells, 
and  d=  1,  separated  by  one  intervening  pixel.  only  vertical  and 
horizontal  separations  were  used. 


The  textural  features  were  calculated  over  each  IB  with  a  differ¬ 
ence  of  1  and  2.  In  addition  the  textural  features  were  calcu¬ 
lated  over  sub  blocks  within  the  IB.  The  means  and  standard 
deviations  of  the  sub  block  statistics  were  saved  for  the  analysis. 
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6  Feature  Discrimination 


The  feature  discrimination  function  is  required  to  separate  class 
types  reliably.  Our  current  classification  function  is  the  Maha- 
lanobis  distance,  which  is  the  multidimensional  analogue  of  the 
Gaussian  distribution  function.  The  category  means  and  crowwcor-, 
relations  are  those  determined  by  Garand  (1986)  and  Garand  and 
We inman ( 1986 ) . 

An  example  of  the  feature  classification  is  the  classification  of 
the  West  Coast  on  Oct  5,  1989.  The  general  region  chosen  for 
c lass i f f icat i on  was  between  30N  and  40  N  in  latitude  and  122W  ard 
128W  in  longitude.  Near  simultaneous  GOES  imagery  of  the  region 
are  shown  for  the  visible  (Figure  2)  and  the  infrared  (Fig. 3). 


The  region  chosen  was  the  area  behind  a  frontal  passage,  charac¬ 
terized  by  a  variety  of  stratus  and  stratocumulus  features. 


Figure  2  -  GOES  Visible  Figure  3  -  IR  image  of  analysis 


Image  of  Analysis 


The  physical  classifier  was  used  to  classify  24  regions  within 
the  image.  Each  IB  was  separated  into  4  32  x  32  pixel  squares 
which  were  used  for  individual  analysis.  The  initial  predictions 
were  found  to  be  highly  skewed  to  clear  classifications.  The 
classifier  was  changed  to  eliminate  clear  values.  The  resulting 
first,  second,  and  third  choices  are  shown  in  Table  4. 

The  visible  image  was  also  analyzed  for  textural  features.  Each 
64  x  64  pixel  region  was  analyzed  for  sum,  difference  and  abso¬ 
lute  difference  vectors  using  the  formulation  of  Eq  2.  The  three 
grey  level  distributions  for  each  distance  were  used  to  find  the 
scene  texture  character istics  of  each  image. 

The  64x64  pixel  array  was  further  divided  into  separate  8X8  pixel 
blocks.  Similar  gray  shade  distributions  were  calculated  for 
each  subregion.  The  mean  and  standard  deviation  of  the  texture 
cha rac te r i s t  ics  were  calculated,  and  are  shown  in  Table  5. 


Table  4 

Physical  Classification  of  DMSP  image 


Ut  N 

Lon  V 

Choice  1 

Choice  2 

Choice  3 

40 .0 

134.9 

Scattered  St 

Dense  CiSt 

Bkn/Ovcst  St 

40.0 

132.0 

Scattered  St 

Broken  Cu 

Dense  CiSt 

jo  .  9 

132.0 

Scattered  St 

Dense  CiSt 

Bkn/Ove*t  St 

49.0 

130.0 

Scattered  St 

Broken  Cu 

Dense  CiSt 

38.0 

130.0 

Scatcered  St 

Broken  Cu 

Dense  CiSt 

38.9 

139.9 

Scattered  St 

Dense  CiSt 

Bkn/Ovcst  St 

34.0 

130.0 

Dense  CISt 

Poly  Open  Ce 

Scattered  St 

48.0 

128.0 

Scattered  St 

Broken  Cu 

Dense  CiSt 

38.0 

128.0 

Scattered  St 

.  Broken  Cu 

Dense  CiSt 

36.9 

128.0 

Dense  CiSt 

Scattered  St 

Poly  Open  Ce 

34.0 

128.0 

Scattered  St 

Dense  CiSt 

Bkn/Ovcst  St 

38.9 

126.0 

Dense  CiSt 

Stratus 

Scattered  St 

36.9 

126.0 

Scattered  St 

Broken  Cu 

Dense  CiSt 

34.8 

126.8 

Scattered  St 

Broken  Cu 

Bkn/Ovcst  St 

48.8 

124.0 

Broken  Cu 

Scattered  St 

Bkn/Ovcst  St 

38.0 

124.0 

Scattered  St 

Broken  Cu 

Bkn/Ovcst  St 

36.8 

124.0 

Scattered  St 

Broken  Cu 

Bkn/Ovcst  St 

34.9 

124.0 

Scattered  St 

Dense  CiSt 

Bkn/Ovcst  St 

36.0 

122.0 

Broken  Cu 

Scattered  St 

Bkn/Ovcst  St 

34.9 

122.0 

Scattered  St 

Dense  CiSt 

Ovrcst  CuNl 

36.0 

120.9 

Poly  Open  Ce 

Dense  CiSt 

Scattered  St 

34.0 

120.0 

Broken  Cu 

Scattered  St 

Table  5 

Textural  characteristics  of  DHSP  visible  image 


